Biophysical Journal

Most kinesin molecular motors dimerize to move processively and efficiently along microtubules; however, some can maintain processivity even in a monomeric state. Previous studies have suggested that asymmetric potentials between the motor domain and microtubules underlie this motility. In this study, we demonstrate that the kinesin-3 family motor protein KLP-6 can move along microtubules as a monomer upon release of autoinhibition. This motility can be explained by a change in length between the head and tail, rather than by asymmetric potentials. Using mass photometry and single-molecule assays, we confirmed that activated full-length KLP-6 is monomeric both in solution and on microtubules. KLP-6 possesses a microtubule-binding tail domain, and its motor domain does not exhibit biased movement, indicating that the tail domain is crucial for the processive movement of monomeric KLP-6. We developed a mathematical model to explain the unidirectional movement of monomeric KLP-6. Our model concludes that a slight conformational change driven by neck-linker docking in the motor domain enables the monomeric kinesin to move unidirectionally if a second microtubule-binding domain exists. SIGNIFICANCEKinesin molecular motors are designed to move efficiently using two heads. Studying these biological molecular motors provides valuable insights into the mechanisms that generate unidirectional movements amidst intense thermal fluctuations. This study reveals that the monomeric kinesin-3 KLP-6 can move along microtubules through interactions with its tail domain. The proposed Brownian ratchet model explains this movement by considering a change in stalk length caused by neck-linker docking rather than asymmetric potentials. This model suggests that a slight conformational change can achieve robust processive movement of kinesin. These findings have significant implications for understanding Brownian ratchet motors and designing rational artificial molecular motors.

For bacterial mechanosensitive channels acting as turgor-adjusting osmolyte release valves, membrane tension is the primary stimulus driving opening transitions. Because tension is transmitted through the surrounding lipid bilayer, it is possible that the presence or absence of different lipid species may influence the function of these channels. In this work, we characterize the lipid dependence of chromosome-encoded MscS and MscL in E. coli strains with genetically altered lipid composition. We use two previously generated strains that lack one or two major lipid species (PE, PG, or CL) and engineer a third strain that is highly enriched in CL due to the presence of hyperactive cardiolipin synthase ClsA. We characterize the functional behavior of these channels using patch-clamp and quantify the relative tension midpoints, closing rates, inactivation depth, and the rate of recovery back to the closed state. We also measure the osmotic survival of lipid-deficient strains, which characterizes the functional consequences of lipid-mediated channel function at the cell level. We find that the opening and closing behavior of MscS and MscL tolerate the absence of specific lipid species remarkably well. The lack of cardiolipin (CL), however, reduces the active MscS population relative to MscL and decreases the closing rate, slightly increasing the propensity of MscS toward inactivation and slowing the recovery process. The data points to the robustness of the osmolyte release system and the importance of cardiolipin for the adaptive behavior of MscS.

Experimental studies reveal that anionic lipid POPA and non-phospholipid cholesterol inhibit the gating of voltage-sensitive potassium (Kv) channels at 5-10% molar concentrations. Intriguingly, other anionic lipids similar to POPA, like POPG, have minimal impact on the gating of the same channels for reasons that remain obscure. Our long-timescale atomistic simulations show that POPA preferentially solvates the voltage sensor domains of Kv channels by direct electrostatic interactions between the positively charged arginine and negatively charged phosphate groups. Cholesterol solvates the voltage sensor domains through CH-{pi} interactions between the cholesterol rings and the aromatic side chains of phenylalanine and tyrosine residues. A continuum electromechanical model predicts that POPA lipids may restrict the vertical motion of voltage-sensor domain through direct electrostatic interactions, while cholesterol may oppose the radial motion of the pore domain of the channel by increasing the mechanical rigidity of the membrane. The electromechanical model predictions are consistent with measurements of the activation curves of Kv channels for various lipids. The atomistic simulations also suggest that the solvation due to POPG is much weaker likely due to its bigger head-group size. Thus the channel activity appears to be tied to the local lipid environment, allowing lipids to regulate channel gating in low concentrations.

During clathrin-mediated endocytosis, a patch of flat plasma membrane is internalized to form a vesicle. In mammalian cells, how the clathrin coat deforms the membrane into a vesicle remains unclear and two main hypotheses have been debated. The "constant area" hypothesis assumes that clathrin molecules initially form a flat lattice on the membrane and deform the membrane by changing its intrinsic curvature while keeping the coating area constant. The alternative "constant curvature" hypothesis assumes that the intrinsic curvature of the clathrin lattice remains constant during the formation of a vesicle while the surface area it covers increases. Previous experimental studies were unable to unambiguously determine which hypothesis is correct. In this paper, we show that these two hypotheses are only two extreme cases of a continuum spectrum if we account for the free energies associated with clathrin assembly and curvature generation. By tracing the negative gradient of the free energy, we define vesiculation pathways in the phase space of the coating area and the intrinsic curvature of clathrin coat. Our results show that, overall, the differences in measurable membrane morphology between the different models are not as big as expected, and the main differences are most salient at the early stage of endocytosis. Furthermore, the best fitting pathway to experimental data is not compatible with the constant-curvature model and resembles a constant-area-like pathway where the coating area initially expands with minor changes in the intrinsic curvature, later followed by a dramatic increase in the intrinsic curvature and minor change in the coating area. Our results also suggest that experimental measurement of the tip radius and the projected area of the clathrin coat will be the key to distinguish between models.

Distance between fluorescent spots formed by various kinetochore proteins ( Delta) is proposed to reflect the level of intrakinetochore tension (IKT). However, larger-scale changes in the kinetochore architecture may also affect Delta. To test this possibility, we measure Delta in long kinetochores of Indian muntjac (IM) whose shape, size, and orientation are discernable in conventional light microscopy. We find that architecture of IM kinetochores and the value of Delta change minimally when microtubule-mediated forces are suppressed by Taxol. In contrast, large decreases of Delta observed in Taxol-treated human cells coincide with prominent changes in length and shape of the kinetochore. We also find that inner and outer kinetochore proteins intermix within a common spatial compartment instead of forming separate thin layers. These observations, supported by computational modelling, suggest that changes in Delta reflect changes in the kinetochore shape rather than the level of IKT.

Inside cells, molecular motors transport cargoes within the actively fluctuating environment known as the cytoplasm. How fluctuations in the cytoplasm affect motor-driven transport is not fully understood. In this study, we investigated the role of fluctuations for transport along microtubules using C. elegans early embryos, focusing on transport driven by cytoplasmic dynein. An artificial motor-cargo complex showed faster transport in vivo than in vitro, suggesting an in vivo acceleration mechanism. We also found that endogenous early endosome transport by dynein is significantly enhanced by the fluctuations in the cytoplasm, which is attributed to the activity of actomyosin networks. An in vitro force measurement of dynein suggests that the asymmetric force response would play a key role in the acceleration. This study provides insights into a regulatory mechanism of molecular motors within actively fluctuating cytoplasm, potentially utilizing random force originating from fluctuating dynamics in the cytoplasm to increase transport efficiency.

Prokaryotic mechanosensitive (MS) channels have an intimate relationship with membrane lipids. Membrane lipids may influence channel activity by directly interacting with bacterial MS channels or by influencing the global properties of the membrane such as area stretch and bending moduli. Previous work has implicated membrane stiffness as a key determinant of the mechanosensitivity of E. coli (Ec)MscS. Here we systematically tested this hypothesis using patch fluorometry of azolectin liposomes doped with lipids of increasing area stretch moduli. Increasing DOPE content of azolectin liposomes causes a rightward shift in the tension response curve of EcMscS. These rightward shifts are further magnified by the addition of stiffer forms of PE such as the branched chain lipid DPhPE and the fully saturated lipid DSPE. Furthermore, a comparison of the branched chain lipid DPhPC to the stiffer DPhPE showed a rightward shift in the tension response curve in the presence of the stiffer DPhPE. We show that these changes are not due to changes in membrane bending rigidity as the tension threshold of EcMscS in membranes doped with PC18:1 and PC18:3 are the same, despite a two-fold difference in their bending rigidity. We also show that after prolonged pressure application sudden removal of force in softer membranes causes a rebound reactivation of EcMscS and we discuss the relevance of this phenomenon to bacterial osmoregulation. Collectively, our data demonstrate that membrane stiffness is a key determinant of the mechanosensitivity of EcMscS.

G protein-coupled receptor signaling has been posited to occur through either collision coupling or pre-assembled complexes with G protein transducers. To investigate the dynamics of G protein signaling, we introduce fluorescence covariance matrix analysis (FCMA), a novel implementation of fluorescence cumulant analysis applied to spectrally resolved fluorescence images. We labeled the GPCR, G, and G{beta}{gamma} units with distinct fluorescent protein labels and we applied FCMA to measure directly the complex formation during stimulation of dopamine and adrenergic receptors. To determine the prevalence of hetero-oligomers, we compared the GPCR data to those from control samples expressing three fluorescent protein labels with known stoichiometries. Interactions between G and G{beta}{gamma} subunits determined by FCMA were sensitive to stimulation with GPCR ligands. However, GPCR/G protein interactions were too weak to be distinguished from background. These findings support a collision coupling mechanism rather than pre-assembled complexes for the two GPCRs studied.

The release of inorganic phosphate (Pi) from the cross-bridge is a pivotal step in the cross-bridge ATPase cycle leading to force generation. It is well known that Pi release and the force-generating step are reversible, thus increase of [Pi] decreases isometric force by product inhibition and increases the rate constant kTR of mechanically-induced force redevelopment due to the reversible redistribution of cross-bridges among non-force-generating and force-generating states. The experiments on cardiac myofibrils from guinea pig presented here show that increasing [Pi] increases kTR almost reciprocally to force, i.e., kTR {approx} 1/force. To elucidate which cross-bridge models can explain the reciprocal kTR-force relation, simulations were performed for models varying in sequence and kinetics of 1) the Pi release-rebinding equilibrium, 2) the force-generating step and its reversal, and 3) the transitions limiting forward and backward cycling of cross-bridges between non-force-generating and force-generating states. Models consisting of fast reversible force generation before/after rapid Pi release-rebinding fail to describe the kTR-force relation observed in experiments. Models consistent with the experimental kTR-force relation have in common that Pi binding and/or force-reversal are/is intrinsically slow, i.e., either Pi binding or force-reversal or both limit backward cycling of cross-bridges from force-generating to non-force-generating states. STATEMENT OF SIGNIFICANCEPrevious mechanical studies on muscle fibers, myofibrils and myosin interacting with actin revealed that force production associated to phosphate release from myosins active site presents a reversible process in the cross-bridge cycle. The correlation of this reversible process to the process(es) limiting kinetics of backward cycling from force-generating to non-force-generating states remained unclear. Experimental data of cardiac myofibrils and model simulations show that the combined effects of [Pi] on force and the rate constant of force redevelopment (kTR) are inconsistent with fast reversible force generation before/after rapid Pi release-rebinding. The minimum requirement in sequential models for successfully describing the experimentally observed nearly reciprocal change of force and kTR is that either the Pi binding or the force-reversal step limit backward cycling.

A broad range of cellular functions involve transient or persistent changes in the morphology of lipid membranes, from the organellar to the molecular scale. By and large, the thermodynamics of these remodeling processes remain to be understood. Molecular Dynamics simulations enhanced by advanced sampling methods are uniquely suited to examine and quantitate these phenomena. Here, we focus on the cellular process known as mechanosensation and use the Multi-Map simulation method to quantify how applied lateral tension impacts the energetics of both global and localized membrane perturbations induced extrinsically. We also examine how tension impacts the dynamics of lipid molecules. We find that the conformational energetics of the membrane clearly differs when it is stretched, and that this difference increases with the magnitude of the applied tension. The reason is not that tension alters the mechanical properties of the lipid bilayer, such as its bending modulus, but rather that it opposes any reduction in the projected area of the membrane relative to that at rest, while the opposite is favored. It follows that tension may shift a conformational equilibrium of a protein that deforms the membrane differently in alternative functional states, if that difference also entails a change in the projected membrane area. Conversely, we find that stretch has little to no effect on the dynamics of lipids at the single-molecule level, implying it would also have no bearing on the lifetime of specific protein-lipid interactions. Finally, we show how changes in lipid composition that result in global membrane thinning can mimic the effect of lateral stretch without any applied tension. Statement of SignificanceCells have evolved the ability to sense mechanical forces, such as pressure or stretch, through specialized proteins embedded in their membranes. How exactly the membrane transduces these stimuli to the proteins therein has been unclear. Using state-of-the-art computer simulations, we show that stretching a membrane does not result in forces that pull or push on the individual lipid molecules that constitute the membrane. Instead, lateral tension alters the energetics of reshaping the membrane. This shift in plasticity explains why several well-known force-sensing proteins switch between active and inactive states at specific tension values observed experimentally. We also show that altering the lipid composition of the membrane can produce the same effect as lateral stretch, without any applied force.

Divalent cations in a concentration-dependent manner behave as effective crosslinkers of intermediate filaments (IFs) such as vimentin IF (VIF). These interactions have been mostly attributed to their multivalency. However, ion-protein interactions often depend on the ion species, and these effects have not been widely studied in IFs. Here we investigate the effects of two biologically important divalent cations, Zn2+ and Ca2+, on VIF network structure and mechanics in vitro. We find that the network structure is unperturbed at micromolar Zn2+ concentrations, but strong bundle formation is observed at a concentration of 100 M. Microrheological measurements show that network stiffness increases with cation concentration. However, bundling of filaments softens the network. This trend also holds for VIF networks formed in the presence of Ca2+, but remarkably, a concentration of Ca2+ that is two orders higher is needed to achieve the same effect as with Zn2+, which suggests the importance of salt-protein interactions as described by the Hofmeister effect. Furthermore, we find evidence of competitive binding between the two divalent ion species. Hence, specific interactions between VIFs and divalent cations are likely to be an important mechanism by which cells can control their cytoplasmic mechanics. SignificanceIntermediate filaments are key structural elements within cells; they are known to form networks that can be crosslinked by divalent cations, but the interactions between the ions and the filaments are not well understood. By measuring the effects that two divalent cations, zinc and calcium, have on the structure and mechanics of vimentin intermediate filaments (VIFs), we show that although both have concentration-dependent effects on VIFs, much more calcium is needed to achieve the same effect as a small amount of zinc. Furthermore, when mixtures of the ions are present, the results suggest that there is binding competition. Thus, cells may use the presence of different cation species to precisely control their internal mechanical properties.

The increase in enzyme-catalyzed reaction rates with temperature is typically modeled using Arrhenius or Eyring relations. Interpretation of extracted parameters is subject to multiple caveats. Here we analyze the impact of expected small temperature variations of underlying activation/Eyring parameters on this modeling. Linear Arrhenius/Eyring behavior can still be observed when the underlying activation energy Ea or enthalpy {Delta}H and entropy {Delta}S vary with temperature. Modest variations -- of the order of an H-bond energy over 60 {degrees}C -- lead to large fractional deviations of Ea, {Delta}H and {Delta}S values derived from linear fits from their underlying values and to deviations of Arrhenius prefactors A by orders of magnitude. In a family of related enzymes with similar activation free energies {Delta}G, small differences in temperature variation of {Delta}H and {Delta}S will lead to apparent enthalpy-entropy compensation and may scramble enzyme ordering based on {Delta}H or {Delta}S. For enzymes from cold and warm-adapted species having largely similar active sites, small temperature variations of {Delta}H and {Delta}S may explain large differences in apparent {Delta}H values. Similar considerations apply to interpretation of van t Hoff plots of equilibrium measurements and related observations of enthalpy-entropy compensation. Complementary methods including simulations and multi-temperature static and time-resolved atomic-resolution structural studies should play a key role in interpreting temperature-dependent kinetic and equilibrium data from enzymatic systems. One-Sentence SummaryLinear Arrhenius/Eyring behavior can occur even when enzyme structure and interactions vary with temperature, yielding fit parameters far from underlying values and generating apparent enthalpy-entropy compensation.

Agonists are evaluated by a concentration-response curve (CRC), with a midpoint (EC50) that indicates potency, a high-concentration asymptote that indicates efficacy and a low-concentration asymptote that indicates constitutive activity. A third agonist attribute, efficiency ({eta}), is the fraction of binding energy that is applied to the conformational change that activates the receptor. We show that {eta} can be calculated from EC50 and the asymptotes of a CRC derived from either single-channel or whole-cell responses. For 20 agonists of skeletal muscle nicotinic receptors, the distribution of {eta} values is bimodal with population means at 51% (including ACh, nornicotine and DMPP) and 40% (including epibatidine, varenicline and cytisine). The value of {eta} is related inversely to the size of the agonists head-group, with high-versus low-efficiency ligands having an average volume of 70 [A]3 versus 102 [A]3. Most binding site mutations have only a small effect on ACh efficiency except for Y190A (35%), W149A (60%) and those at G153 (42%). If {eta} is known, the midpoint and high-concentration asymptote can be calculated from each other. Hence, an entire CRC can be estimated from the response to a single agonist concentration, and efficacy can be estimated from EC50 of a CRC that has been normalized to 1. Given {eta}, the level of constitutive activity can be estimated from a single CRC. Statement of significanceReceptors are molecular machines that convert chemical energy from agonist binding into mechanical energy of a global conformational change that generates a cell response. Agonists are distinguished by their potency (proportional to affinity) and efficacy, but also by the efficiency at which their binding energy is applied to receptor activation. Here we show that agonist efficiency can be estimated from a single concentrationresponse curve, and estimate efficiencies of 20 nicotinic receptor agonists. These have a bimodal distribution with larger agonists belonging to the lower-efficiency population. We further show that mutations of some binding site residues alter efficiency, and that knowledge of agonist efficiency simplifies and extends dose-response curve analysis.

Outer hair cells (OHCs) are essential for the sensitivity and frequency specificity of the mammalian ear. To perform this function, OHCs need to amplify the motion of the basilar membrane, which is much stiffer than themselves. OHCs must overcome this impedance mismatch for their amplifying function particularly at high frequencies, where the mismatch is largest. This issue could be solved by the existence of multiple modes of motion. Here, systems of two coupled oscillators are examined as the simplest of such cases. It is found that some of these model systems indeed make OHCs function as an effective amplifier by overcoming the impedance mismatch. This result suggests that the elaborate structure of the organ of Corti, which can support multiple modes of motion, is a key to the high frequency performance of the mammalian ear. SignificanceThe mammalian ear depends on outer hair cells, which work as the cochlear amplifier. The mechanism, with which outer hair cells perform this biological function, is of great interest. The present paper addresses a question, how soft outer hair cells can amplify the vibration of the much stiffer basilar membrane. It shows that the elaborate structure of the cochlea, which supports multiple modes of motion, must be a key to the exquisite performance of the mammalian ear. It also shows that the properties of outer hair cells obtained from isolated cell preparations are compatible with their physiological function.

Studying the role of molecularly distinct lipid species in cell signaling remains challenging due to a scarcity of methods for performing quantitative lipid biochemistry in living cells. We have recently used lipid uncaging to quantify lipid-protein affinities and rates of lipid transbilayer movement and turnover in the diacylglycerol signaling pathway using population average time series data. So far, this approach does not allow to account for the cell-to-cell variability of cellular signaling responses. We here report a framework that allows to uniquely identify model parameters such diacylglycerol-protein affinities and transbilayer movement rates at the single cell level for a broad variety of structurally different diacylglycerol species. We find that lipid unsaturation degree and longer side chains generally correlate with faster lipid transbilayer movement and turnover and higher lipid-protein affinities. In summary, our work demonstrates how rate parameters and lipid-protein affinities can be quantified from single cell signaling trajectories with sufficient sensitivity to resolve the subtle kinetic differences caused by the chemical diversity of cellular signaling lipid pools.

The accumulation of viral structural proteins along the ER-Golgi intermediate compartment (ERGIC) membrane leads to SARS-CoV-2 self-assembly and budding, driven by the interactions between these proteins, RNA and the ERGIC membrane. The membrane protein (M) is believed to interact with other structural proteins and form clusters needed for the induction of membrane curvature that facilitates virion formation. However, the role played by direct and membrane-mediated interactions between M proteins and their interactions with other proteins in the clustering process remains unclear. Here, we utilize a combination of all-atom molecular dynamics (MD) simulations, continuum modeling and experiments to show that M-M interactions are sufficient to drive clustering in ERGIC-like lipid bilayers in the absence of other proteins or RNA. Using all-atom MD simulations we were able to estimate the membrane thinning induced by M proteins and the resulting membrane-mediated M-M interaction. Combining this with a continuum model that describes the evolution of M protein density in a planar lipid membrane, we identified the existence of a critical, direct M-M interaction energy needed for cluster assembly at a given density. By comparing the model predictions with analysis of atomic force microscopy images of M protein clusters in supported lipid bilayers, we were able to estimate the direct M-M interaction energy and found it to be significantly larger than the membrane mediated interaction energy. Our work therefore establishes that M protein interactions are sufficient to drive clustering and provides a quantitative understanding of the role played by direct and membrane-mediated interactions of M proteins in viral assembly and budding.

Model asymmetric bilayers are useful for studying the coupling between lateral and transverse lipid organization. Here, we used calcium-induced hemifusion to create asymmetric giant unilamellar vesicles (aGUVs) for exploring the phase behavior of 16:0-PC/16:1-PC/Cholesterol, a simplified model for the mammalian plasma membrane. Symmetric GUVs (sGUVs) were first prepared using a composition that produced coexisting liquid-disordered and liquid-ordered phases visible by confocal fluorescence microscopy. The sGUVs were then hemifused to a supported lipid bilayer (SLB) composed of uniformly mixed 16:1-PC/Cholesterol. The extent of outer leaflet exchange was quantified in aGUVs in two ways: (1) from the reduction in fluorescence intensity of a lipid probe initially in the sGUV ("probe exit"); or (2) from the gain in intensity of a probe initially in the SLB ("probe entry"). These measurements revealed a large variability in the extent of outer leaflet exchange in aGUVs within a given preparation, and two populations with respect to their phase behavior: a subset of vesicles that remained phase separated, and a second subset that appeared uniformly mixed. Moreover, a correlation between phase behavior and extent of asymmetry was observed, with more strongly asymmetric vesicles having a greater probability of being uniformly mixed. We also observed substantial overlap between these populations, an indication that the uncertainty in measured exchange fraction is high. We developed models to determine the position of the phase boundary (i.e., the fraction of outer leaflet exchange above which domain formation is suppressed) and found that the phase boundaries determined separately from probe-entry and probe-exit data are in good agreement. Our models also provide improved estimates of the compositional uncertainty of individual aGUVs. We discuss several potential sources of uncertainty in the determination of lipid exchange from fluorescence measurements. Statement of SignificanceWe used calcium-induced hemifusion to create an asymmetric lipid distribution in giant unilamellar vesicles that are models for the mammalian plasma membrane. Confocal fluorescence micrographs of asymmetric vesicles showed that coexisting liquid-ordered and liquid-disordered domains initially present in symmetric vesicles were disrupted after 75% of the saturated lipid in their outer leaflets was replaced with unsaturated lipid. We developed quantitative models for extracting valuable information from the data, including the location of the phase boundary and the compositional uncertainty of individual asymmetric vesicles. The methodology we describe can help reveal the molecular determinants of interleaflet coupling of phase behavior and thus contribute to a better understanding of lipid raft phenomena.

The Transient Receptor Vanilloid 1 (TRPV1) is a non-selective ion channel, which is activated by several chemical ligands and heat. We have previously shown that activation of TRPV1 by different ligands result in single-channel openings with different conductance, suggesting that the selectivity filter is highly dynamic. TRPV1 is weakly voltage-dependent, here we sought to explore whether the permeation of different monovalent ions could influence the voltage-dependence of this ion channel. By using single-channel recordings, we show that TRPV1 channels undergo rapid transitions to closed states that are directly connected to the open state, which may result from structural fluctuations of their selectivity filters. Moreover, we demonstrate that the rates of these transitions are strongly influenced by the permeant ion, suggesting that ion permeation regulates the voltage dependence of these channels.

We address the challenge of understanding how hydrophobic interactions are encoded by fusion peptide sequences within coronavirus (CoV) spike proteins. Within the fusion peptides of SARS-CoV-2 and MERS-CoV, a largely conserved peptide sequence called FP1 (SFIEDLLFNK and SAIEDLLFDK in SARS-2 and MERS, respectively) has been proposed to play a key role in encoding hydrophobic interactions that drive viral-host cell membrane fusion. While a non-polar triad (LLF) is common to both FP1 sequences, and thought to dominate the encoding of hydrophobic interactions, FP1 from SARS and MERS differ in two residues (Phe 2 versus Ala 2 and Asn 9 versus Asp 9, respectively). Here we explore if single molecule force measurements can quantify hydrophobic interactions encoded by FP1 sequences, and then ask if sequence variations between FP1 from SARS and MERS lead to significant differences in hydrophobic interactions. We find that both SARS-2 and MERS wild-type FP1 generate measurable hydrophobic interactions at the single molecule level, but that SARS-2 FP1 encodes a substantially stronger hydrophobic interaction than its MERS counterpart (1.91 {+/-} 0.03 nN versus 0.68 {+/-} 0.03 nN, respectively). By performing force measurements with FP1 sequences with single amino acid substitutions, we determine that a single residue mutation (Phe 2 versus Ala 2) causes the almost threefold difference in the hydrophobic interaction strength generated by the FP1 of SARS-2 versus MERS, despite the presence of LLF in both sequences. Infrared spectroscopy and circular dichroism measurements support the proposal that the outsized influence of Phe 2 versus Ala 2 on the hydrophobic interaction arises from variation in the secondary structure adopted by FP1. Overall, these insights reveal how single residue diversity in viral fusion peptides, including FP1 of SARS-CoV-2 and MERS-CoV, can lead to substantial changes in intermolecular interactions proposed to play a key role in viral fusion, and hint at strategies for regulating hydrophobic interactions of peptides in a range of contexts. SIGNIFICANCEFusion of coronaviruses (CoVs) and host cells is mediated by the insertion of the fusion peptide (FP) of the viral spike protein into the host cell membrane. Hydrophobic interactions between FPs with their host cell membranes regulate the viral membrane fusion process and are key to determining infection ability. However, it is not fully understood how the amino acid sequences in FPs mediate hydrophobic interactions. We use single-molecule force measurements to characterize hydrophobic interactions of FPs from SARS-CoV-2 and MERS-CoV. Our findings provide insight into the mechanisms by which the amino acid composition of FPs encodes hydrophobic interactions and their implications for fusion activity critical to the spread of infection.

Molecular motors convert chemical potential energy into mechanical work and perform a great number of critical biological functions. Examples include the polymerization and manipulation of nucleic acids, the generation of cellular motility and contractility, the formation and maintenance of cell shape, and the transport of materials within cells. The mechanisms underlying these molecular machines are routinely divided into two categories: Brownian ratchet and power stroke. While a ratchet uses chemical energy to bias thermally activated motion, a stroke depends on a direct coupling between chemical events and motion. However, the multi-dimensional nature of protein energy landscapes allows for the possibility of multiple reaction paths connecting two states. Here, we investigate the properties of a hypothetical molecular motor able to utilize parallel ratchet and stroke translocation mechanisms. We explore motor velocity and force-dependence as a function of the energy landscape of each path and reveal the potential for such a mechanism to result in an optimum force for motor function. We explore how the presence of this optimum depends on the rates of the individual paths and show that the distribution of stepping times characterized by the randomness parameter may be used to test for parallel path mechanisms. Lastly, we caution that experimental data consisting solely of measurements of velocity as a function of ATP concentration and force cannot be used to eliminate the possibility of such a parallel path mechanism. SIGNIFICANCEMolecular motors perform various mechanical functions in cells allowing them to move, replicate and perform various housekeeping functions required for life. Biophysical studies often aim to determine the molecular mechanism by which these motors convert chemical energy to mechanical work by fitting experimental data with kinetic models that fall into one of two classes: Brownian ratchets or power strokes. However, nothing a priori requires that a motor function via a single mechanism. Here, we consider a theoretical construct where a motor has access to both class of mechanism in parallel. Combining stochastic simulations and analytical solutions we describe unique signatures of such a mechanism that could be observed experimentally. We also show that absence of these signatures does not formally eliminate the existence of such a parallel mechanism. These findings expand our theoretical understanding of the potential motor behaviors with which to interpret experimental results.